Direct-space investigation of the ultraslow ballistic dynamics of a soft glass 
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We use light microscopy to investigate the aging dynamics of a glass made of closely packed soft 
spheres, following a rapid transition from a fluid to a solid-like state. By measuring time-resolved, 
coarse-grained displacements fields, we identify two classes of dynamical events, corresponding to 
reversible and irreversible rearrangements, respectively. The reversible events are due to the small, 
experimentally unavoidable fluctuations of the temperature imposed to the sample, leading to tran- 
sient thermal expansions and contractions that cause shear deformations. The irreversible events 
are plastic rearrangements, induced by the repeated shear cycles. We show that the displacement 
due to the irreversible rearrangements grows linearly with time, both on average and at a local 
level. The velocity associated with this ballistic motion decreases exponentially with sample age, 
accounting for the observed slowing down of the dynamics. The displacement field due to the irre- 
versible rearrangements has a vortex-like structure and is spatially correlated over surprisingly long 
distances. 

PACS numbers: 82.70.-y, 61.20.Lc, 61.43.-j, 62.20.Fe 



I. INTRODUCTION 

A large variety of soft disordered materials, ranging 
from colloidal pastes and gels, hard-spheres suspensions, 
dense surfactant, copolymer or star polymer systems, ex- 
hibit very slow dynamics, often associated with physical 
aging, the continuous growth of relaxation times after 
the sample has been prepared In this respect, these 
systems display intriguing analogies with molecular or 
polymer glasses and are often termed "Soft Glassy 
Materials" Analogies are also found with granular 
materials UH, B0] close to the jamming [1] transition: 
these similarities have spurred a large research activity, 
aiming at rationalizing the glass and jamming transitions 
within one single conceptual framework Q. Note how- 
ever that the driving force for the slow dynamics is dif- 
ferent in molecular and colloidal systems, where ther- 
mal energy dominates, and in athermal granular media, 
where some form of external drive is required to induce 
structural rearrangements. 

Among soft glasses, a class of systems have emerged 
whose dynamics exhibit peculiar features not seen nei- 
ther in molecula r g lass formers nor in granular me- 
dia 0, M M [H, M M M M [13, [H, m • They range 
from tenuous, diluted gels made of attractive particles, 
to repulsive colloidal glasses or concentrated surfactant 
phases. In these systems, correlation functions measured 
by dynamic light scattering or X-photon correlation spec- 
troscopy decay steeper than exponentially ( "compressed 
exponential" ) and their decay time scales as the inverse 
scattering vector, q, rather than as g"^, as expected for a 
diffusive process. This suggests that the displacement of 
the scatterers grows linearly with time, e.g. that the slow 
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dynamics is, on average, ballistic. Recent time-resolved 
light scattering experiments [l^l on a tenuous colloidal 
gel, the first system for which these unusual dynamics 
were reported [9| , have shown that particle displacements 
are not continuous but rather result from intermittent re- 
arrangements that affect a very large region of the sam- 
ple [2l|. 

In spite of the large number of scattering experiments 
reporting ultraslow ballistic relaxations, the ballistic na- 
ture of the dynamics and its discontinuous behavior have 
never been demonstrated by directly i mag ing the motion 
of particles. Direct space information [22| on the dynam- 
ics would be very valuable in order to better understand 
the physical origin of such unusual dynamics. It is gen- 
erally proposed that ballistic dynamics stem from the re- 
laxation of internal stress HI] . While a plausible mi- 
croscopic picture of the mechanisms through which stress 
may be accumulated and released has been proposed (but 
not directly demonstrated) for tenuous colloidal gels, a 
convincing explanation for concentrated systems is still 
lacking. In particular, it is unclear whether thermal en- 
ergy alone could be responsible for significant restructur- 
ing in very dense colloidal or surfactant systems, where 
the elastic modulus can be a factor of lO'^ — 10"' larger 
than that of the tenuous gels. 

To address some of these issues, we have recently used 
light microscopy to investigate the slow dynamics of a soft 
glass made of a compact arrangement of elastic spheres 
("onions") [l^], the same system for which ballistic dy- 
namics were reported in previous light scattering experi- 
ments [ll| . In Ref. [131 we used a coarse-grained method 
to measure two-time displacement fields in direct space. 
The main conclusion of this work was that the slow dy- 
namics is due to shear deformations resulting from the 
experimentally unavoidable small fluctuations of the tem- 
perature imposed to the sample. Indeed, as temperature 
fluctuates the sample undergoes a series of contractions 
and elongations, because of thermal expansion. While 
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macroscopically these deformations are fully reversible, 
microscopically they eventually induce irreversible rear- 
rangements. These experiments suggested that the slow 
dynamics of dense colloidal systems may share more simi- 
larities than expected with sheared athermal systems. In 
Ref. [ij we mainly focussed on average quantities, such 
as the mean squared displacement and its age depen- 
dence. Here, we extend this study and provide a novel 
and more direct method to quantify the irreversible re- 
arrangements, showing unambiguously that they lead to 
ballistic motion, both on average and at a local level. 
Furthermore, we investigate in detail the spatial struc- 
ture of the irreversible rearrangements. We find that the 
displacement field is correlated over distances comparable 
to the system size, presenting interesting analogies with 
recent simulations [25i. [261] of jammed spheres and exper- 
iments on very concentrated colloidal suspensions [27[ . 

The paper is organized as follows. Section HIl describes 
the experimental system, the measurement procedure 
and the image analysis method that we have developed. 
In Sect ion Hm we first present the time-resolved (but spa- 
tially averaged) dynamics and then discuss the individual 
trajectories and the spatial structure of the rearrange- 
ments. In the last section, we discuss briefly our main 
findings and conclude. 



II. MATERIALS AND METHODS 
A. Sample preparation and optical microscopy 

The structure and the phase diagram of our system 
have been discussed in detail in Ref. [28| : here we simply 
recall their main features. The sample is a surfactant 
lamellar phase constituted of a regular stacking of bilay- 
ers. The bilayers are composed of a mixture of cetylpyri- 
dinium chloride (CpCl) and octanol (Oct) (weight ratio 
CpCl/Oct = 0.95) and diluted in brine ([NaCl] = 0.2M) 
at a weight fraction (p = 16%. An amphiphilic copolymer, 
Symperonics F68 by Serva ((EO)76 - (PO)29 - (EO)76, 
where EO is ethylene oxide and PO is propylene oxide) is 
incorporated into the lamellar phase. At room tempera- 
ture the PO block is hydrophobic and adsorbs onto the 
bilayers, while the hydrophilic EO blocks remain swollen 
in the aqueous solvent and decorate the surfactant bi- 
layers. Upon copolymer addition, we have shown both 
experimentally and theoretically that a transition from a 
flat lamellar phase to an onion phase occurs [1^ . The lat- 
ter consists of a dense packing of multilamellar vesicles, 
or onions. Due to their polydispcrsity and softness, the 
volume fraction of the onions is one. In the experiments 
reported here, the copolymcr-to-bilaycr weight ratio is 
0.8 and the radius of the largest onions is about 6 /im 

The rheological properties of our samples can be tuned 
by changing the temperature, T. Both the EO and PO 
blocks of the F68 copolymer are hydrophilic at low tem- 
perature (T < 8°C): the copolymer is fully soluble in 



water and the sample is fluid. Because the PO block be- 
comes increasingly hydrophobic when increasing T, the 
sample becomes gel-like at room temperature, where the 
PO block adsorbs onto the surfactant lamellae. The 
temperature behavior provides a convenient and well- 
reproducible way to initialize the dynamics, by rapidly 
quenching the sample from its low-T fluid phase to the 
solid phase pU. [23.l29|. At room temperature, the onion 
phase behaves mechanically as a gel, with a storage mod- 
ulus, G' ~ 300 Pa, that is nearly frequency-independent 
in the range 5 x 10^"^ — 10 Hz and roughly one order of 
magnitude larger than the loss modulus, G" . 

We probe the dynamics of the onion phase by means 
of optical microscopy. The sample is contained in a 
glass capillary of rectangular cross-section with thickness 
200 /im and width 2 mm, which is flame-sealed in order 
to prevent evaporation. Centrifugation is used to confine 
the air bubble at one end of the capillary (see Fig. [T|). 
The sample is placed in an oven (Instec) that controls 
the temperature T to within 0.15°C (standard deviation 
of T over 1 day). T is measured with a thermocouple 
placed in direct contact with the capillary. A light mi- 
croscope equipped with a lOx objective is used to image 
a thin slice of the sample, in the mid plane between the 
upper and lower wall of the capillary. The field of view 
is 0.93 mm x 1.24 mm and is located in the center of the 
capillary, as shown in the scheme of Fig. [1] The onions 
are difficult to visualize using conventional optical mi- 
croscopy, due to their weak refractive index mismatch 
with respect to the solvent, and because they are densely 
packed. While individual onions in diluted samples are 
clearly visible using differential interference contrast mi- 
croscopy [13, closely-packed onions are best visualized 
under crossed polarizers. Due to their birefringence and 
their spherical shape, they appear as dark crosses on a 
bright background, as can be seen in Fig. [TJ We use a 
charge-coupled detector camera (CCD) to take images of 
the sample, with a delay between images of 15 sec; each 
image results from an average over 10 frames, in order to 
reduce the electronic and read-out noise of the CCD. 

The sample is loaded in the capillary in its fluid state 
and the capillary is immediately placed under the micro- 
scope, in the oven regulated at 5°C. The dynamics is 
then initialized by a T jump from 5°C to 23.3°C, which 
is performed in situ. The time needed to reach the fi- 
nal temperature is less than one minute. We define age 
tu, = as the time at which T has reached its stationary 
value 23.3°C. A typical experiment lasts about 1 day: 
we have checked that no measurable drift of the capillary 
due to mechanical instabilities occurs over this time. 



B. Image Analysis 

Due to the tight packing and to polydispcrsity, a typ- 
ical onion image consists of a highly contrasted inten- 
sity pattern, where only a few, very large onions may 
be individually resolved, as shown in Fig. [TJ Thus, it is 




FIG. 1: (Top) Scheme of the capillary, with the air bubble, 
the location of the field of view (dotted area), and the ori- 
entation of the parallel and perpendicular axis. The size of 
the field of view is 0.93 mm x 1.24 mm. (Bottom) Portion of 
a typical image of the sample as taken by light microscopy 
between cross-polarizers. The arrow points to a large onion: 
note the characteristic four-lobed pattern due to the onion 
birefringence. 

impossible, in general, to track the motion of individual 
objects. Instead, we use a method inspired by the parti- 
cle imaging velocimetry (PIV) technique widely used in 
fluid mechanics to obtain coarse-grained displacements 
fields of fluids seeded with tracer particles Simi- 
lar approaches specialized to the case where individual 
particles are not discernable, as in our case, have been 
termed digital imaging correlation (DIG) [31| or image 
correlation velocimetry (ICV) js^]. In these methods, a 
time series of images of the sample is taken, usually using 
a CCD camera. Each image is divided into a grid of re- 
gions of interest (ROIs) , corresponding to square regions 
of side L in the sample. The internal dynamics within a 
given ROI is supposed to be negligible compared to the 
overall drift of that ROI. In other terms, the variation 
of the displacement field over a distance L is supposed 
to be much smaller than its absolute magnitude. Under 
this assumption, the intensity pattern of a given ROI in 
an image taken at time t will appear to be shifted in the 
corresponding ROI imaged at time t + t. The displace- 
ment can be estimated by calculating by what amount a 
ROI of the second image has to be back-shifted in order 
to maximize its resemblance with the corresponding ROI 
of the first image. By repeating this procedure for all 
ROIs, a coarse-grained displacement field is obtained. 

Usually, for a given ROI the shift along the direction 
of columns {5x) and rows ((5y), is determined by search- 
ing the maximum of the spatial crosscorrelation of the 
intensity of the two images. With this method, the dis- 
placement is determined with a resolution of one pixel. 
Several schemes have been proposed to improve this res- 
olution, e.g. by calculating the center of mass of the peak 



of crosscorrelation, or by fitting it to a 2-dimensional an- 
alytical function such as a Gaussian. We find that these 
standard approaches do not perform very well in our case, 
most likely due to the anisotropic shape of the correla- 
tion peak resulting from the peculiar intensity pattern in 
the onion images. We have therefore developed an alter- 
native approach based on a least-square method, where 
we minimize x^{5x, 5y ), the variance of the difference be- 
tween the first image, /, and a shifted version of the sec- 
ond image, J. In order to calculate the optimum shift 
with subpixel resolution, we use linear interpolation and 
define as follows: 

X\5x,5y) =j:r,cKciSx,Sy) (1) 

ArASx, Sy) =Ir^c-{l-5x)il-Sy)Jr^c- (2) 
Sx{l - Sy)Jr,c+i - (1 - Sx)SyJr+i,c 
~Sx5yJr+i,c+i ■ 

Here, Ir^c is the intensity at time t of the pixel at row 
r and column c and J^.c the corresponding quantity at 
time t + T. The sum is over all pixels belonging to a given 
ROI. For simplicity, Eq. (2) has been written for the case 
< Sx < 1, < Sy < 1. The general case and the details 
of the algorithm will be discussed in a forthcoming paper. 

We have performed several tests of our algorithm. 
First, we have checked that the gradient of the displace- 
ment field is not too high, as required by PIV-based 
methods. We find that the difference of the displacement 
field measured for adjacent ROIs is typically at most 20% 
of the displacement itself, thus validating our approach. 
Furthermore, we find that different choices of the ROI 
size yield consistent results, as long as L is at least a 
few onion diameters. As one would expect, increasing 
the ROI size reduces the measurement noise, but results 
in a coarser displacement field. The best trade-off be- 
tween high spatial resolution and low noise is found for 
L ~ 78 yum, corresponding to 47 pixels. The full dis- 
placement field is then composed of iV = 192 ROIs, on 
a 12 rows x 16 columns grid. To quantify the accuracy 
and the precision of our measurements, we have used a 
piezoelectric motion controller to repeatedly translate a 
sample by a known amount (1 /im). The measured dis- 
placement field was found to agree with the imposed one 
to within 80 nm in the worst case. In the following, we 
will take this value as the typical uncertainty of our mea- 
surements. 



III. RESULTS 

A. Time-resolved dynamics 

The raw data obtained from the image analysis are 
two-time displacement maps, /S.Wj{t^,T). Here, a 
stands either for || or _L and refers to the component 
of the displacement along or perpendicular to the cap- 
illary long axis (see Fig. [1]). The superscript i,j desig- 
nate the i-th row and the j-th column of the ROI grid 
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(1 < i < 12 and I < j < 16), and the displacements are 
measured between images taken at time and + r. 
Although in the following we will discuss the behavior of 
both components, the focus will be on Ai?y'-' (iw, t), since 
the dynamics along the long axis of the capillary typically 
dominates the one in the perpendicular direction. 

In addition to the local dynamics, we also analyze time- 
resolved global quantities, obtained by averaging over 
space. For each component, we define ARa{tw,T), the 
spatially averaged absolute displacement between time 
tw and tw + t: 



Ai?„(i„,T) = l^Ai?j;^(t„,T), 

id 



(3) 



where the sum is over the N ROIs of the grid. We also 
introduce the corresponding spatially averaged relative 
displacement: 




(4) 



where Ar^J(i„,T) = Ai?^^(t„,T) - Ai?„(t„,T) is the 
local deviation with respect to the spatially averaged dy- 
namics. The relative displacements are by definition pos- 
itive and quantify the spatially heterogeneous character 
of the instantaneous displacement field. The quantities 
defined above are function of both the age of the sample, 
tji,, and the lag between images, r. In Ref. [l^l Ar|| (tto, t) 
has been studied as a function of t^' here, we fix the age 
of the sample and explore the evolution of the absolute 
and relative displacements as a function of the lag r be- 
tween images. Although most results presented in the 
following refer to a specific age, we have checked that 
they are representative of the general behavior of the 
sample. 

We plot in Fig. [2^ the lag dependence of the parallel 
absolute displacement, Ai?|| (i^,, r), for an age = 15390 
s and a restricted, representative range of lags. Ai?|| is 
temporally heterogenous, with intermittent peaks, whose 
magnitude is of the order of 1 /im. The time between 
these events is of the order of a few hundreds seconds. 
Neither the amplitude of the peaks nor their frequency 
appear to evolve systematically with r. Indeed, the 
displacement averaged over a time window significantly 
longer than the typical time between peaks, but signifi- 
cantly shorter than the duration of the experiment, does 
not vary with time (data not shown). This indicates 
that the sample does not drift, and that the dynamics 
is stationary. In Ref. HJ] we have shown that the in- 
termittent behavior is due to small, experimentally un- 
avoidable fluctuations of the temperature imposed by the 
oven. These fluctuations arc of the order of 0.1°C and 
lead to elongations and contractions of the sample, due 
to the thermal expansions of water, the main constituent 
of our system. The fluctuations of the spatially aver- 
aged absolute displacement Ai?|| are quite peculiar. They 
exhibit abrupt negative-going swings, corresponding to 
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FIG. 2; (Color on line) The three panels show, for a sample 
age tu, = 15 390 s and as a function of lag r, the spatially- 
averaged absolute displacement parallel to the capillary axis 
(a) and the spatially-averaged relative displacement parallel 
and perpendicular to the capillary axis (b and c, respectively). 
The black lines are the data for all lags 500 > r < 7400 s, 
while the cyan triangles indicate irreversible events alone, as 
defined in the text. 



AT < 0, followed by smoother positive displacements. 
These features reflect the time evolution of T due to the 
oven thermostat, which cools abruptly while heating up 
more smoothly. Because of the capillary geometry, T 
fluctuations lead to an uniaxial elongation or contraction 
of the sample, along the || axis. In agreement with a uni- 
axial thermally- induced expansion/compression mecha- 
nism, we measure that the absolute, spatially averaged, 
displacement along the _L axis is always very small and 
comparable to the measurement uncertainty (data not 
shown) . 

We show in Fig. [^b and [2t the parallel and perpendic- 
ular spatially averaged relative displacements, Ay'n and 
Ar±, as a function of the lag, for the same lag window 
as that used in Fig. The small range of lags shown 
in Fig. [5] allow a close comparison of the three signals. 
The perpendicular relative displacement (Fig. O;) in- 
creases smoothly with time, without any significant tem- 
poral fluctuation. By contrast, the parallel component of 
the relative displacement (Fig. ^p) exhibits intermittent 
peaks, whose amplitude is several tenths of /im. This in- 
termittent behavior is similar to that of the absolute par- 
allel displacement shown in a). In [23| we have analyzed 
the same quantities for a fixed lag, as a function of t^j, 
and have shown that the intermittent peaks of the rela- 
tive and absolute displacement are simultaneous and that 
both are induced by temperature fluctuations. The same 
conclusions arc reached here; note that both positive and 
negative peaks of AR arc associated to positive spikes of 
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Ar, since the latter is by definition positive. The analysis 
of AR and Ar vs r presented here further clarifies the 
nature of the dynamics of the onions: while the absolute 
displacement fluctuates around a constant value [sst , the 
parallel and perpendicular relative displacements overall 
increase with lag t. In other words, even when the sam- 
ple comes back to a previous position on average, as in- 
dicated by AR = 0, the local configuration has changed, 
since Ar ^ 0: this proves directly the existence of irre- 
versible rearrangement events. It is worth noting that 
the scenario emerging from these observation is remi- 
niscent of the dynamics of dense emulsions or colloidal 
suspensions to which a periodic shear is applied [13, UHl , 
although the physical origin of the motion is different. In 
those systems, the particle positions decorrelate as a re- 
sult of an applied shear stress. In our sample, elongations 
and contractions due to temperature fluctuations arc not 
spatially uniform, as we will discuss in more detail in 
Sec. nil Bl Therefore, they result in a shear deformation 
analogous to that experienced by the systems in [sl, HI] . 
In the experiments described in [3^ . l35j , not all emulsion 
droplets or colloidal particles recover their initial position 
when the applied stress is removed, because some plastic 
rearrangements have occurred. Similarly, the configura- 
tion of our system is permanently changed after several 
contraction/elongation cycles, because irreversible rear- 
rangements have occurred. 

In order to characterize the irreversible rearrangement 
events, for a given we identify the lags for which the 
absolute parallel displacement is null. For these lags, we 
expect no contribution to the displacement field due to 
the shear induced by contraction or elongation, since the 
sample is globally undeformed. Hence, any local displace- 
ment can be attributed exclusively to irreversible rear- 
rangements. In practice, we identify the lags t for which 
|Ai?||| < (0 ± 0.08) /im, where the range of accepted ab- 
solute displacements is chosen equal to the measurement 
uncertainty. These lags are marked by triangles in Fig. [2l 
We find that all the symbols lie on the "baseline" of Ar|| , 
indicating that this procedure indeed identifies only those 
lags for which the change in sample configuration is min- 
imal, i.e. lags for which no extra contribution due to a 
transient shear is present. 

We show in Fig. [3] the lag-dependence of the relative 
displacements both parallel, Arj"', and perpendicular, 

Ar^'' , to the long axis of the capillary for the irreversible 
events alone. Different symbols correspond to data anal- 
ysis performed for four close but distinct sample ages in 
the range 15030 s < tw < 16140 s. The four sets of 
data nicely collapse onto two unique curves. Remark- 
ably, we find that, up to r w 35000 s, i.e. over more 
than a decade in time, Arj" varies linearly with the lag. 

A fit with a straight line, Ar]"' = Vb(tw)''' + b, is shown 
in Fig. [Si. We find b = (0.13 ± 0.1)nm, of the order of 
the measurement uncertainty. Thus, up to experimental 
uncertainties, the initial growth of Ary is proportional to 
t: the ultraslow motion associated with irreversible rear- 
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FIG. 3: (Color on line) r dependence of the spatially-averaged 
relative displacements, parallel (empty symbols) and perpen- 
dicular (full symbols) to the capillary long axis, for the irre- 
versible events alone, as defined in the text. The four different 
symbols correspond to data taken at close but different sam- 
ple ages: tw = 15030 s (squares), 15390 s (circles) , 15780 s 
(up triangles) and 16140 s (down triangles). Lin-lin plot in a) 
and log-log plot in b). In b) an offset on the order of 0.13 /im 
has been subtracted to the data. The black line in a) is a 
linear fit of the data for r < 35000 s. In b), the black line is 
a power law with an exponent of 1. 



rangements is ballistic. For the data of Fig. [31 taken at 
sample ages from — 15030 s to = 16140 s, we find 
Vh = (3.8±0.1) X lO-^ms-^ The ballistic nature of the 
dynamics for not too long lags is robust: indeed, ballistic 
motion is measured for all sample ages for lags shorter 
than T « 35000 s. As shown in Fig. 31 Vb decreases 
decreases steadily with sample age, An exponential 
decay, with a characteristic decay time of (72000 ± 3000) 
s, can fit the data, although experimental measurements 
over a larger range of sample ages would be required to 
better characterize the functional form of the decrease. 
In Ref. [2^ the contribution of the irreversible events to 
the overall dynamics has been estimated using a different 
method, which led to similar conclusions: the irreversible 
dynamics is ballistic, with a characteristic velocity that 
decays nearly exponentially with sample age. The de- 
cay time of Vb found here is about twice as large as the 
one (40000 ± 5500 s) evaluated in [U, probably due to 
the difference in the analysis method. The analysis of 
AR and Ar vs r presented above shows that the typi- 
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times tw + Ti and + ti + T2 (see inset of Fig. [61d) . More 
specifically, we calculate 
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FIG. 4: (Symbols) Ballistic velocity for irreversible motion 
along the parallel direction as a function of sample age. The 
error bars are the standard deviations for analysis performed 
for at least 4 close but distinct sample ages. The line is an 
exponential fit yielding a decay time of (72000 ± 3000) s. 



cal frequency of the contracions/elongations causing the 
irreversible rearrangements does not change with i^, as 
expected for a mechanism driven by externally imposed 
temperature fluctuations. We thus conclude that the ag- 
ing revealed by the decrease of 14 must be due to a de- 
crease of the amplitude of the displacement associated 
with each event. 

The departure of Arj" from a straight line for lags 
larger than k, 35000 s (see Fig. ^ is due to the slowing 
down of the dynamics, and the decrease of the velocity 
with sample age. The slowing down of the dynamics is 
more pronounced in the perpendicular direction than in 
the parallel direction. Indeed, a saturation of Ar'['' at 
a value of about 0.6 /im is observed for lags larger than 
~ 35000 s, while no saturation is observed for ArJ" over 
the duration of the experiment. 



B. Individual trajectories and spatial structure of 
the dynamics 

The ballistic motion discussed so far has been inferred 
from the behavior of l^rit^j^r), a quantity that is spa- 
tially averaged. To better characterize the irreversible 
events, we analyze also the time evolution of the local 
relative displacements. To this end, we study the trajec- 
tory of individual ROIs. where only delays corresponding 
to irreversible events are considered. Representative tra- 
jectories located at different positions in the grid of ROIs 
are shown in Fig. [51 where the displacements relative to 
the position at = 15390 s are plotted as color-coded 
circles, according to the lag r. We find that all trajecto- 
ries are close to straight lines, the main changes between 
different ROIs being the length of the trajectory and its 
orientation, i.e. the magnitude and the orientation of the 
local velocity. 

To quantify our observations, we evaluate for all ROIs 
the correlation between the direction of the relative dis- 
placement between times and iw+Ti and that between 
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Ar'-J(tw,n) ■ Ar'-J(tw + ri,T2) 
|Ar'J(iw,Ti)||Ar*.J'(tw +Ti,r2)| 



(5) 



and its average over all ROIs, cos6'. If the trajectories 
were random walks, cos 9^'^ would be distributed between 
— 1 and 1, and its average value would be close to zero. By 
contrast, for ballistic motion we expect cos 9 close to 1. 
We plot in Fig.[6K the value of cos 9, in the plane (r2,ri). 
The color-coded map shows that, except for small ti or 
small T2 for which results are dominated by experimental 
uncertainties, the value of cos6' is always larger than 0.5, 
suggesting a strong correlation of the direction of motion 
over long time scales. 

We now consider the case of two successive displace- 
ments with approximately the same duration, hence T2 ~ 
Ti . We thus show in Fig. we show cos 9 for T2 ~ n . 
The data are obtained from the points in Fig. [5^ that lie 
in between the two straight lines T2 = ri ± 20% shown 
in the figure. At short lags, for which Arj" and Ar^'' 
are not significantly different from the experimental un- 
certainties, cos9 varies randomly between —0.5 and 0.5 
before reaching rapidly a plateau of height 0.67 ±0.06 for 
Ti > 4500 s, corresponding to an average angle of 48 dcg. 
The maximum value for cos 6* is « 0.8, for ri ~ 10000 
s (corresponding to an angle 9 of 37 deg). The slightly 
decreasing trend above 10000 s is due to the saturation 
of Ar^"^ at large r (Fig. [3]), which is associated with a 
kink in individual trajectories, as shown in Fig. O Over- 
all, cos 9 is close to unity, indicating that the direction of 
motion persists over long times (as long as the duration 
of the experiment). This analysis demonstrates that the 
motion is ballistic not only when averaged over the whole 
field of view, but also at a local level, since most individ- 
ual ROIs move in a well defined direction and exhibit 
trajectories close to straight line over the whole duration 
of the experiment. 

What is the spatial structure of the displacement field 
associated with such ballistic dynamics? By examining 
a large number of displacement fields for all ages and all 
delays, we find that generally they can be divided in two 
classes. The first one corresponds roughly to a shear de- 
formation along the long axis of the capillary, as shown in 
Fig. [7Ji.. This displacement pattern is associated with a 
significant elongation or contraction of the sample; thus 
it corresponds to the transient, reversible peaks of the 
relative displacement. The second class is characterized 
by a vortex-likc (or compression-like) structure, as seen 
in Fig.[7}3. Note that, since the sample is actually incom- 
pressible, Fig.lTjj suggests that some motion may occur in 
the direction perpendicular to the plane that is imaged. 
In contrast to the shear-like deformation, this pattern is 
quite isotropic and implies non-negligible displacements 
along the _L direction. The events that have been pre- 
viously identified as irreversible are typically associated 
with a relative displacement field of the second kind. The 
overall pattern of the displacement fields associated with 
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FIG. 5: (Color online) Trajectory of 4 single ROIs with re- 
spect to the position at time = 15390 s. The location of 
each ROI on the 12 x 16 overall grid is specified by the row 
and column indices i and j. Only data points correspond- 
ing to irreversible events as defined in the text are plotted. 
The color of the circles changes according to the delay r, as 
shown by the scale. The spatial scale is the same for the four 
trajectories. 



irreversible events does not significantly vary with time, 
as expected from the analysis of the individual trajecto- 
ries, which are all very close to straight lines. 

From our analysis, we can conclude that the T- 
induccd, intermittent events and the ballistic irreversible 
rearrangements can be associated to distinct, specific 
structures of the relative displacement field. The former 
correspond to the (anisotropic) shear deformation of the 
sample along the capillary axis (Fig. Uh.), while the lat- 
ter correspond to nearly isotropic, compression-like dis- 
placement fields (Fig.lTb). Both structures are correlated 
over distances comparable to the field of view. We have 
checked that they do not exhibit significant changes with 
sample age, as expected from the ballistic motions of the 
individual ROIs. 



IV. DISCUSSION AND CONCLUSION 

We have used optical microscopy to study the slow 
dynamics of closely packed onions, a soft glassy mate- 
rial. Two classes of dynamical events have been identi- 
fied: reversible and irreversible rearrangement. We have 
shown that the reversible, intermittent peaks of relative 
displacement are due to the contraction or elongation of 
the sample, due to temperature fluctuations, and corre- 
spond to a shear deformation along the long axis of the 
capillary. The origin of the shear associated with ther- 
mal expansion or contraction remains somehow unclear. 
It may be due to the curvature of the meniscus at the 
interface between the free end of the sample and the air 
bubble. It may also stem from small thermal gradients 
in the sample, close to the capillary wall. Indeed, por- 
tions of the sample close to the lateral wall are in better 
thermal contact with the oven, with respect to the cen- 
tral region. Therefore, they may react faster to any small 
change in T and lead to a larger elongation or contrac- 
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tion, as we observe in our experiments. Wc thus expect 
the precise spatial structure of the displacement fields to 
depend on the sample geometry (e.g. to be different in 
a microscopy capillary, a light scattering cuvette, and a 
Couettc cell for rhcology measurements), but the effect 
of T fluctuations should be relevant in all experimental 
configurations. 

The second class of dynamical events are the irre- 
versible rearrangements, presumably occurring as the re- 
sult of repeated shear cycles. They correspond to a nearly 
isotropic displacement field that is spatially correlated 
over very large distances. The motion resulting from 
these ultraslow rearrangements is ballistic, with a veloc- 
ity that decreases exponentially as the sample ages. We 
believe that these irreversible events occur at all times, 
presumably even when they are masked by larger, tran- 
sient shear events caused by temperature fiuctuations. 

Below, we discuss shortly the main features of the irre- 
versible rearrangements and finally comment on the role 
of the temperature fluctuations as the driving force for 
the slow dynamics. 
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FIG. 7: a) Representative example of a shear-like relative 
displacement field, observed for r = 615 s and tu, = 15390 
s. b) Vortex-like relative displacement field, observed for r = 
68100 s and t„ = 15390 s. The field of view is 0.93 mm x 
1.24 mm; the arrows are magnified 50 times. 



A. Ballistic motion and aging 

We have shown that the amplitude of the relative dis- 
placements is proportional to the lag, r, hence prov- 
ing directly the ballistic nature of the ultraslow motion. 
A similar ballistic motion has been invoked to explain 
the slow dynamics of a variety of soft glassy materi- 
als pro bed by light and X-photon scattering techniques 
(Tll.1T^. [13I . [14 . [la . [isl . [36j including the very same onion 



sample investigated here [ll[ . To the best of our knowl- 



edge, our experiments are the first proof in direct space 
of ballistic motion in the slow relaxation of a glassy ma- 
terial. Below, we provide simple scaling arguments to 
understand the order of magnitude of the ballistic veloc- 
ity measured experimentally. 

In Ref. [2§|, we have presented a model to rationalize 
the dynamical behavior of the onion glass, speculating 
that motion results from the combined action of viscous 
forces and internal elastic stress. We briefly recall the 
main features of the model: regions of size D containing 
several onions move under the action of elastic forces due 
to the relaxation of internal stresses and experience a 
viscous drag. The driving force due to the relaxation of 
internal stress is = crl?^, where a = eG is the elastic 
stress, G the elastic modulus, and e the strain. Stokes' 



law gives the viscous force: F^j = pr]Dv where p is a 
geometrical factor of order 10 that depends on the shape 
of the region, rj the effective viscosity of the medium, and 
V the velocity of the region. By balancing the viscous and 
the elastic force and solving for w, one obtains 



V = GeD/iiip) . 



(6) 



We calculate v using experimental values for 77, e, G and 
D. To estimate the order of magnitude of the strain, 
we consider the irreversible displacement fields measured 
experimentally and take the difference of the relative dis- 
placement between a ROI for which it is close to and 
one for which it is maximum (about 1 fim for i„ = 15000 
s) . The typical distance between these two ROIs is about 
500 ^m, yielding e « 0.002. In our microscopy measure- 
ments, we probe the motion of ROIs of size L, so we set 
D = L ~ 78 //m. We take for 77 the value of the effective 
viscosity deduced from step strain relaxation measure- 
ments [23], 77 = 3.4 X 10^ Pas for an age t„ = 15000 s, 
and take G ~ 300 Pa as measured by linear rheology [1^ . 
Using Eq. we then find v « 10^^"/p ms^^ with p 
typically of the order of ten. This result is in very good 
agreement with Vb{t^ = 15000 s) = 3.8 x 10"" ms'^ 
measured by light microscopy. 

It is interesting to compare Vb to Vols, the ballistic 
velocity obtained in [2^ from the analysis of DLS data 
for the same sample. For an age — 15000 s we found 
in [ll| Vdls = 5 X 10"^'^ ms~^, about two orders of mag- 
nitude lower than Vt, measured, at the same age, by light 
microscopy. Several factors may be invoked to explain 
the large difference. First, the dynamics is not probed on 
the same length scale in the two experiments: here, the 
relevant length scale is L = 78 /xm, the ROI size. In the 
DLS experiment, by contrast, the motion is probed on a 
much smaller length scale, set by the inverse scattering 
vector and of the order of 0.2 — 1 /xm. Furthermore, the 
microscopy experiments point to the role of T-induced 
strain fields as the driving force for the slow rearrange- 
ments. In this scenario, the latter most likely vary with 
both the frequency and the amplitude of the transient 
strain, which depend on experimental conditions such as 
the cell size and geometry and the temperature control, 
different in microscopy and DLS measurements. More 
experiments, e.g. varying the cell geometry, would be 
useful to clarify this issue. 

The slowing down of the dynamics with sample age, 
a feature typical of glassy dynamics, is reflected here by 
the decrease of Vb with t„. 14 is a global quantity, av- 
eraged over both space and time. Its decay may a priori 
stem from one or more distinct mechanisms: a decrease 
of the rate at which rearrangements occur, a reduced av- 
erage displacement associated with each rearrangement, 
or a smaller size of the volume affected by a single rear- 
rangement. Our space- and time resolved measurements 
rule out the former and the latter hypothesis: rearrange- 
ments occur at the same rate, imposed by the characteris- 
tic time of the temperature fluctuations of the oven, and 
affect essentially the same (large) portion of the sam- 
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pie at all ages. Experiments on concentrated colloidal 
hard spheres [37| have reached a somehow similar con- 
clusion: the size of the "clusters" formed by the fastest 
moving particles did not evolve significantly during aging 
(Note that these clusters contained just a few particles, 
in contrast to the very long-ranged dynamical correla- 
tions observed here). A few simulation studies have ad- 
dressed the issue of the age dependence of the correlation 
length of the dynamics, yielding a different answer: the 
correlation length was found to grow with both in a 
Lennard- Jones [38l. [soj and a spin glass [13, El . More ex- 
periments and simulations will be needed to understand 
these discrepancies. 



B. Long-ranged spatial correlations of the slow 
dynamics 

We have shown that the relative displacement fields as- 
sociated with the irreversible events exhibit surprisingly 
long-ranged spatial correlations, comparable to the size 
of the field of view, thus much larger than the onions' 
size. This result is strikingly different from previous ex- 
perimental and numerical determinations of the size of 
dynamical heterogeneity, which generally are of the or- 
der of a few particle sizes, e.g. in dense colloidal hard 
spheres [H, El, granular materials 0, EE Ell, and 
molecular glass formers [H, E3 El, Ell- Recent sim- 
ulations in a two-dimensional hard sphere system may 
however reconcile these findings The authors 

show that the rearrangements are short-range correlated 
for supercooled, equilibrated samples, while they be- 
come correlated over long distances deep in the out-of- 
equilibrium glass phase, approaching the maximum pack- 
ing. While previous works where short-ranged correla- 
tions were found generally dealt with equilibrium super- 
cooled samples, it is clear that our closely-packed, aging 
onions are jammed systems, for which long -ranged cor- 
relations may be expected according to f25r[26j . Indeed, 
we have recently observed very long-ranged spatial cor- 
relations in the dynamics of a variety of jammed systems 
using space-resolved dynamic light scattering [2l|, con- 
sistently with the scenario emerging here. 

As a further intriguing analogy, we note that the dis- 
placement field following a rearrangement event in p6j 
has a vortex-likc structure that is reminiscent of that 
observed in our experiments. Similar structures have 
been also reported for the non-affine part of the displace- 
ment measured, experimentally and numerically, in vari- 
ous amorphous systems under shear, ranging from granu- 



lar material to glasses [50|, |51l 153, l53|, |5j] . Whether these 
analogies are accidental or not remains an open question. 



C. Role of the temperature fluctuations 

The experiments described in this paper highlight the 
key role of temperature fiuctuations in determining the 
slow dynamics of the onions. One may wonder whether 
T-induced shear deformations are the only driving force 
for the non-stationary irreversible rearrangements. In- 
deed, it is difficult to imagine that thermal energy alone 
could promote such a large scale motion. In fact, the 
onions are a closely-packed, fully jammed system with a 
solid-like behavior. In this respect, they are similar to 
dry foams or compressed emulsions, none of witch ex- 
hibits a spontaneous slow dynamics in the absence of a 
(non-thermal) driving force, such as an applied stress or 
the stress continuously build up by coarsening (ssl . |56| . 
For jammed or nearly jammed systems, on the other 
hand, an applied shear strain or stress has been shown to 
promote irreversible rearrangements in a variety of sys - 
tems, including granular media 0, H, emulsions [s^ . 
and concentrated suspensions of colloidal [s^l and non- 
colloidal [H, [ill particles. 

A final question concerns the generality of the mech- 
anism reported here. In our time-resolved light scat- 
tering experiments on dilute colloidal gels we did 
not observe any correlation between temperature fluctu- 
ations and the rearrangement events that were detected. 
By contrast, a strong correlation has been established 
for the onion glass and has also been observed in time- 
resolved light scattering measurements on a concentrated 
micellar phase 
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We speculate that temperature- 
induced strain fluctuations may be important in very con- 
centrated, solid-like systems where expansions and con- 
tractions necessarily displace all constituents, and where 
thermal energy alone may be insufficient to significantly 
rearrange the structure. On the other hand, in dilute 
systems such as the gels, expansions and contractions of 
the solvent would be accommodated fiowing through the 
porous structure, without affecting the particles. More 
experiments with a finer temperature control or an exter- 
nally imposed shear will be necessary to test these ideas. 

We thank M. Wyart for useful discussions. This work 
was supported in part by the European NoE "SoftComp" 
(NMP3-CT-2004-502235) and the French ACI JC2076. 
L.C. is a junior member of the Institut Universitaire de 
France, whose support is gratefully acknowledged. 



[1] L. Cipelletti, L. Ramos, J. Phys.: Condens. Matter 17 
R253 (2005). 

[2] E. Donth, The Glass Transition (Springer, Berlin, 2001). 
[3] P. Sollich, F. Lequeux, P. Hebraud, M. E. Gates , Phys. 
Rev. Lett. 78, 2020 (1997). 



[4] O. Pouliquen, M. Belzons, and M. Nicolas, Phys. Rev. 

Lett. 91, 014301 (2003). 
[5] G. Marty and O. Dauchot, Phys. Rev. Lett. 94, 015701 

(2005). 

[6] P. M. Reis, R. A. Ingale, and M. D. Shattuck, Phys. Rev. 



10 



Lett. 98, 188301 (2007). 
[7] A. S. Keys, A. R. Abate, S. C. Glotzer, D. J. Durian, 
Nat. Phys. 3, 260 (2007). 

A. J. Liu and S. R. Nagel, Nature 396, 21 (1998). 
L. Cipelletti, S. Mauley, R. C. Ball, and D. A. Weitz, 
Phys. Rev. Lett. 84, 2275 (2000). 

A. Knaebel, M. Bellor, J.-P. Munch, V. Viasnoff, F. 
Lequeux, and J. L. Harden, Europhys. Lett. 52, 73 [34 

(2000) . 

L. Ramos and L. Cipelletti, Phys. Rev. Lett. 87, 245503 [35 

(2001) . 

L. Cipelletti, L. Ramos, S. Mauley, E. Pitard, D. A. [36 
Weitz, E. E. Pashkovski, and M. Johansson, Faraday Dis- 
cuss. 123, 237 (2003). [37 
R. Baudyopadhyay, D. Liang, H. Yardimci, D. A. Ses- 
soms, M. A. Borthwick, S. G. J. Mochrie, J. L. Harden, [38 
and R. L. Leheny, Phys. Rev. Lett. 93, 228302 (2004). [39 

A. Robert, E. Waudersmau, E. Dubois, V. Dupuis, and 
R. Perzynski, Europhys. Lett. 75, 764 (2006). [40 
R. Baudyopadhyay, D. Liaug, J. L. Harden, and R. L. 
Leheuy, Solid State Commuu. 139, 589 (2006). [41 

B. Chung, S. Ramakrishuan, R. Baudyopadhyay, D. 
Liaug, C. F. Zukoski, J. L. Harden, and R. L. Leheuy, [42 
Phys. Rev. Lett. 96, 2006). [43 
P. Falus, S. Narayanan, A. R. Sandy, and S. G. J. 
Mochrie, Phys. Rev. Lett. 97, 066102 (2006). [44 
H. Guo, J. N. Wilkiug, D. Liaug, T. G. Mason, J. L. 
Harden, and R. L. Leheny, Phys. Rev. E 75, 041401 [45 
(2007). 

C. Caronua, Y. Chushkin, A. Madseu, and A. Cupaue, [46 
Phys. Rev. Lett. 1, 055702 (2008). 

A. Duri and L. Cipelletti, Europhys. Lett. 76, 972 (2006). [47 
A. Duri, PhD thesis, Uuiversite Montpellier 2 (2006). [ 
V. Prasad, D. Semwogerere, and E. R. Weeks, J. Phys.: [49 
Coudeus. Matter 19, 113102 (2007). 
J.-P. Bouchaud and E. Pitard, Eur. Phys. J. E 6, 231 
(2001). [50 
S. Mazoyer, L. Cipelletti, and L. Ramos, Phys. Rev. Lett. 
97, 238301 (2006). [51 
C. Brito and M. Wyart, Europhys. Lett. 76, 149 (2007). 
C. Brito and M. Wyart, J. Stat. Mech. L08003 (2007). [52 
P. Ballesta, A. Duri, and L. Cipelletti, Nature Physics 4, 
550 (2008). [53 
L. Ramos, D. Roux, P. D. Olmsted, and M. E. Gates, [54 
Europhys. Lett. 66, 888 (2004). 

L. Ramos and L. Cipelletti, Phys. Rev. Lett. 94, 158301 [55 
(2005). 

G. E. Willert and M. Gharib, Exp. Fluids 10, 181 (1991). [56 
J. Westerweel, Meas. Sci. Technol. 8, 1379 (1997). 
P. T. Tokumaru and P. E. Dimotakis, Exp. Fluids 19, 1 [57 
(1995). 

[33] AR compares the sample configurations at time tw and [58 
tw + T. If the first image considered (at t^) has been 
taken at high temperature (with respect to the average [59 
temperature), then the reference sample configuration is 
expanded, and AR will be on average more compressed 
than the reference state and will thus fluctuate with time 



around a negative value. By contrast if the first image 
considered has been taken at relatively low temperature, 
then the reference sample configuration is "compressed" , 
and AT? will be on average more expanded than the ref- 
erence state and will thus fluctuate with time around a 
negative value. If the reference state corresponds to a 
relaxed state, then AR will fluctuate around 0. 
P. Hebraud, F. Lequeux, J. P. Munch, and D. J. Pine, 
Phys. Rev. Lett. 78, 4657 (1997). 

G. Petekidis, A. Moussai'd, and P. N. Pusey, Phys. Rev. 
E 66, 051402 (2002). 

M. Bellour, A. Knaebel, J. L. Harden, F. Lequeux, and 
J.-P. Munch, Phys. Rev. E 67, 031405 (2003). 
R. E. Courtland and E. R. Weeks, J. Phys.: Coudeus. 
Matter 15, S359 (2003). 

H. E. Castillo and A. Parsaeian, Nat. Phys. 3, 26 (2007). 
A. Parsaeian, H. E. Castillo, preprint 
|arXiv:c ond-mat /0610789 

J. Kisker, L. Sauten, M. Schreckenberg, and H. Rieger, 
Physical Review B 53, 6418 (1996). 

L. Berthier and J. P. Bouchaud, Physical Review B 66, 
054404 (2002). 

L. Berthier et ai, Science 310, 1797 (2005). 

E. R. Weeks, J. C. Crocker, A. C. Levitt, A. Schofleld, 

D. A. Weitz, Science 287, 627 (2000). 

E. R. Weeks, J. C. Crocker, D. A. Weitz, J. 
Phys.:Condens. Matter 19, 205131 (2007). 

O. Dauchot, G. Marty, G. Biroh, Phys. Rev. Lett. 95, 
265701 (2005). 

F. Lechenault, O. Dauchot, G. Biroli, and J. P. 
Bouchaud, Europhys. Lett. 4, 46003 (2008). 

M. D. Ediger, Annu. Rev. Phys. Chem. 51, 99 (2000). 

S. C. Glotzer, J. Nou-Cryst. Solids 274, 342 (2000). 

C. Dalle- Ferrier, C. Thibierge, C. Alba-Simiouesco, L. 

Berthier, G. Biroli, J.-P. Bouchaud, F. Ladieu, D. LHote, 

and G. Tarjus, Phys. Rev. E 76, 041510 (2007). 

A. Tanguy, J. P. Wittmer, F. Leouforte, and J.-L. Barrat, 

Phys. Rev. B 66, 174205 (2002). 

A. Tanguy, F. Leouforte, and J.-L. Barrat, Eur. Phys. J. 
E 20, 355 (2006). 

E. Kolb, J. Cviklinski, J. Lauuza, P. Claudiu, and E. 
Clement, Phys. Rev. E 69, 031306 (2004). 

C. E. Maloney, Phys. Rev. Lett. 97, 035503 (2006). 

E. Del Gado, P. Ilg, M. Kroger, and H. C. Ottinger, Phys. 
Rev. Lett. 101, 095501 (2008). 

D. J. Durian, D. A. Weitz and D. J. Pine, Science 252, 
686 (1991). 

S. Coheu-Addad and R. Hohler, Phys. Rev. Lett. 86, 
4700 (2001). 

V. Viasnoff and F. Lequeux, Phys. Rev. Lett. 89, 065701 
(2002). 

D. J. Pine, J. P. GoUub, J. F. Brady, and A. M. Leshan- 
sky. Nature 438, 997 (2005). 

L. Corte, P. M. Chaikiu, J. P. GoUub, D. J. Pine, J. 

F. Brady, and A. M. Leshausky, Nature Physics 4, 420 
(2008). 



